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Abstract 

The relativistic amplitude for the direct knock-out contribution to 
(7,p) reactions on nuclei is reduced to a nonrelativistic form using an 
effective Pauli reduction scheme. The reduction is carried out to sec¬ 
ond order in the inverse nucleon mass. It is found that the interaction 
Hamiltonian appearing in the nonrelativistic amplitude has significant 
dependence, starting at second order, on the vector and scalar mean 
nuclear potentials. These strong medium modifications are absent 
in traditional nonrelativistic calculations. Detailed comparisons show 
that these modifications are crucial to understanding the differences 
between relativistic and nonrelativistic models. These differences are 
also examined through reduction of the relativistic amplitude via the 
Foldy-Wouthuysen transformation. Similar medium modifications are 
obtained in this case as well. We discuss the implications of these 
medium modifications for the consistency of existing nonrelativistic 
calculations. 

*Work supported in part by the Natural Sciences and Engineering Research Council of 
Canada 
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Introduction 


Recently Lotz and Sherif have calculated the relativistic amplitude for the 
direct knock-out contribution to ( 7 ,p) reactions in a distorted wave approx¬ 
imation |I|. Similar calculations have also been given by McDermott et al. 
[§]. The results obtained from these relativistic calculations were found to 
be in better agreement with experimental data than those of nonrelativis- 
tic distorted wave Born approximation (DWBA) calculations 0. The latter 
calculations are based on the use of Schrodinger wave functions to describe 
the bound and continuum nucleons. The nuclear current is obtained through 
a Foldy-Wouthuysen (FW) transformation of the relativistic Hamiltonian 
describing the interaction of a photon with a free nucleon |^, 

In this paper we clarify the main differences between the relativistic and 
nonrelativistic calculations. We obtain the nonrelativistic amplitude from the 
relativistic distorted wave S matrix through two different reduction schemes. 
The hrst of these is the Pauli reduction in which the relativistic S matrix 
of the ( 7 ,p) reaction is expressed in a form that involves Schrodinger-like 
wave functions and an effective interaction Hamiltonian. This Hamiltonian 
is expanded in powers of 1/ {E + M), where M is the nucleon mass and E is 
its total energy. The hrst order terms in this expansion, in the limit E M, 
reproduce the form of the usual nonrelativistic amplitude. A characteris¬ 
tic feature of the higher order terms is the dependence of the interaction 
Hamiltonian on the nuclear vector and scalar potentials. A short account of 
this discussion has been published 0. This scheme has also been applied to 
a similar study of the relativistic and nonrelativistic descriptions of (e, e'p) 
reactions 0. 

The second approach is based on the use of the Foldy-Wouthuysen trans¬ 
formation. The relativistic knock-out amplitude is written for a model case 
in which the initial bound and hnal continuum nucleons are described by a 
single Dirac Hamiltonian with real vector and scalar potentials. Through the 
FW reduction we are able to write a nonrelativistic limit of the amplitude to 
various orders in the inverse nucleon mass. Again we hnd that the interaction 
Hamiltonian is dependent on the strong nuclear potentials. 

We begin section by summarizing the calculation of the relativistic S 
matrix describing the direct knock-out contribution to ( 7 ,p) reactions on nu¬ 
clei. In addition we provide some discussion of the Dirac equation containing 
strong scalar and vector potentials. In section ^ we illustrate the formalism 
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of the Pauli reduction and show how the relativistic S matrix of the {'y,p) 
reaction is reduced to an expansion in the inverse nucleon mass. We then 
show results of detailed calculations and discuss the implications of the pres¬ 
ence of the nuclear potentials in the interaction Hamiltonian. In section we 
discuss the FW transformation of a relativistic Hamiltonian describing a par¬ 
ticle interacting with electromagnetic and strong nuclear helds, and we hnd 
the corresponding nonrelativistic amplitude to second order in the inverse 
nucleon mass. We discuss differences between the interaction Hamiltonians 
obtained through the Pauli and FW reduction schemes. The section closes 
with a discussion of observables calculated at different orders, with and with¬ 
out potentials in the FW interaction Hamiltonian. Section ^ is devoted to 
conclusions. 


1 Relativistic Direct Knock-out Mechanism 


The relativistic distorted wave amplitude (S matrix) for the reaction {'y,p) on 
a target nucleus in the angular momentum state | JjMj) leading to a residual 
nuclear state \JfMf) is obtained in hrst order in the interaction Hamiltonian 

as m 
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( 1 ) 


where Ec is the energy of the hnal state continuum proton and M is its mass. 
The energy of the incident photon is oj. The Clebsch-Gordan coefficient 
in equation (|) couples the bound proton with angular momentum Jb and 
projection Mb to the residual nucleus, and S is the spectroscopic factor. The 
electromagnetic interaction Hamiltonian is 

Hem (a:) = (4 {x) + ■ (2) 

Notice that in this work we do not consider modihcations to the electromag¬ 
netic interaction Hamiltonian arising when a nucleon is off-shell |@]. The 
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4-vector potential describing the photon is written as 

Al (a:) = (3) 

where the label ^ is the polarization angle of the photon, allowing us to 
choose between two linear polarization states when the polarization 4-vector 
is written as = (0, cos^, 0). In writing the polarization vector this 
way we have made a dehnite choice of z-axis; namely, the z-axis is taken along 
the direction of the photon momentum so the photon 4-momentum is always 
kt!j = 0, 0, k-y). We thus have two linear polarization states: = 0° with 

polarization along the x-axis and ^ = 90° with polarization along the y-axis. 

The second term in the interaction Hamiltonian (2) is an anomalous mag¬ 
netic moment term with k = 1.79 for the proton. The tensor is related 
to the usual Dirac gamma matrices through, 

^ b", 7l , (4) 

and the electromagnetic held tensor is 

- duAy,. (5) 

The Dirac spinors in equation (|I]) describe the bound and continuum 
nucleons and are solutions of a Dirac equation of the form 


{cx - p + I3[M + S (r)] + V {r)}ip (x) = Eijj (x), 


( 6 ) 


where we adopt the standard representation of the 4x4 Dirac matrices {«*} 
and /3 0. The potentials S (r) and V (r) are the scalar and zeroth-component 
vector potentials respectively. These potentials are originally derived from 
the effective Lagrangian of the cx-l-a; model [|T3- In actual distorted wave cal¬ 
culations, the hnal state proton is described by a solution of the Dirac equa¬ 
tion above, containing complex potentials Sc (r) and Vc (r). The parameters 
for these potentials are obtained through analysis of data from proton elastic 
scattering on nuclei |^. For the bound nucleon, Dirac-Hartree potentials 
are used |T^. Detailed discussions of the Dirac equations pertaining to both 


the bound and continuum nucleons can be found in references 10 


In the rest of this section we discuss a Schrodinger-like wave equation 
derived from the Dirac equation (^. The solutions of this Schrodinger-like 
equation will be used in the next section to cast the amplitude in a nonrela- 
tivistic form. 
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We write the Dirac spinors in terms of upper and lower components u 
and £ 


-0 (x) 


'm (x)' 


(7) 


and use the Dirac equation (P) to write the lower component of the wave 
function in terms of the upper component as 


£ (x) 


__ / N 

E + M + S{r)-V{r)'^ 


( 8 ) 


Thus the Dirac spinors can be written in terms of their upper components 
in the form 


'll) (x) 


1 

cr • p 

M + E + S{r) 



(9) 


The upper component of the Dirac spinor u (x) can be related to a 
Schrodinger-like wave function Tsch (x) by []^ 


u{x) = D 2 (r) Tsch (x), 


( 10 ) 


and the function D (r) depends on the Dirac potentials as 


D{r) 


E + M + S{r)-V{r) 
E + M 


s (r) 

E + M' 


( 11 ) 


Note that D{r) goes to one for large r, so the asymptotic behavior of the 
upper component of the Dirac wave function and the Schrodinger-like wave 
function are the same. The Schrodinger-like wave function Tsch {x) is the 
solution of the Schrodinger-like equation [|T^] 


2M 


+ Ucent {r) + f/so (r) cr ■ T ^ ^sch = {E- M) Tgch- 


( 12 ) 


The central and spin-orbit potentials are functions of the Dirac potentials as 
well as the energy of the proton, and are written explicitly as 


Ucent 


E-M + 
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(13) 
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and 


1 s' (r) 
2M rs (r) 


( 14 ) 


In addition to the fnnction s(r) defined in eqnation (pd]), we have defined a 
fnnction involving the difference of the nncleon’s mass and energy 


d{r) = M-E + S{r) + V{r). 


(15) 


For the bonnd state this “nonrelativistic” eqnation provides a better de¬ 
scription of spin orbit splitting than the usnal nonrelativistic calcnlations 
r^ , and similarly for the continnnm nucleon the wave function obtained 


from this equation gives an improved description of nucleon-nucleus elastic 
scattering data 0. 


2 Effective Pauli Reduction 


Pauli reduction of the relativistic amplitude for the knock-out contribution 
to ( 7 ,p) reactions was described briefly in reference |^; we provide more 
details and discussion here. In the effective Pauli reduction scheme the rela¬ 
tivistic distorted wave S matrix is rewritten in terms of Schrodinger-like wave 
functions resulting in an effective interaction Hamiltonian which may be 
expanded in powers of The various orders can then be related, in the 

proper limit, to the nonrelativistic form of the amplitude. As we show below 
there are, however, important differences with the usual nonrelativistic am¬ 
plitude. The qualitative features of a similar two component reduction of the 
Dirac wave function for several generic vertices, in the presence of the nuclear 

HI 


interactions, has been discussed by Cooper et al. 


The emphasize in our 


present discussion is to present a more quantitative analysis for the case of 
( 7 ,p) reactions. 


2.1 Formalism 

Beginning with the amplitude of equation (|l|), we write the integral in terms 
of Schrodinger-like wave functions using equations (||) and (pT|) relating the 
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Dirac spinors to the Schrodinger-like wave functions. This allows us to write 

-i r 1 1 1/2 r n 1/2 

xl 1/2 
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where the effective interaction Hamiltonian is 


( 16 ) 
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This can be expanded in powers of and written in the form; 

= //}!) + /7j2) ^ _ 

The hrst and second order contributions are given by 
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where we have written 


Qx{r) 


Sxjr) - Vx(r) 
E + M 


( 20 ) 


O [/]■■■ in the interaction Hamiltonians of equation ( P!U[ ) means that operator 
O acts only on function /. The appearance of the Dirac potentials through 
Q’s in the interaction Hamiltonians (ED delineates the modihcation of the 
effective photon interaction due to the presence of the nuclear medium. We 
will investigate the signihcance of this medium effect in the following section. 

Using the Hamiltonian ( 0 ) in equation (|T^ ) along with the nonrelativistic 
wave functions, we can cast the amplitude in a nonrelativistic form as follows 


c<(h _ 

^fi - 


(27r 


i r 1 1 r 1 1/2 
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( 21 ) 


where {i) refers to the highest order of the inverse of the nucleon mass in 
the interaction Hamiltonian used in the nonrelativistic amplitude. The am¬ 
plitude (pTD, with H} = in the limit Eb,Ec —M is equivalent to the 
usual nonrelativistic transition amplitude [^3] except that the Schrodinger- 
like wave functions are used instead of the usual nonrelativistic Schrodinger 
wave functions. The amplitude obtained this way will be referred to here as 
the hrst order nonrelativistic amplitude. 

The Schrodinger-like wave functions describing the bound nucleon can be 
written as 


while for the continuum nucleon we write 

4'Lh,o(2^) = 

LMJ 

X (L, 1/2; M - S,, S/l J, M)hj(r)yt[%Aa), (23) 

where 

= i:(i.i/2;M-^.,,|j,M)yy-'‘(!2)x5-/2. ( 24 ) 
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It should be emphasized here that the wave functions introduced in equations 
(p2D and (p3D have nonrelativistic normalization, i.e the factor \J ^ 2 m 'w^hich 
comes from the Dirac spinor describing the outgoing nucleon has been set 
equal to one (thus a bound state wave function is normalized to one and the 
plane wave limit of the nucleon wave function is of the form d'sch,Sj-(^) = 

addition the factor which comes from the Dirac held 
expansion is set equal to one in the cross section. 

Using equations (|), (p^) and (^) then after evaluating the angular in¬ 
tegration the hrst order amplitude can be written as 


0 ( 1 ) _ 

- 


-le 


r 11 1/2 


TT 


.2a; 

x5{Ec- 


JbMb 


SjiJj{JB) 


1 I /2 


ILJ^ 


X (L, 1/2; Mb + 2/i - Sf, Sf | J, Mb + 2/r)y^' 


Ms+^fi—Sf 


(%) 


X 


KU) 

2M 


+ 


00 


M + Ec 

— (2/r coS(^ — i sin,^) ^ 


3 '^/ T 


M + Ec M + E^ 


B 


^ ^LJLb 


{-^Lb+ 1 , 1 ,L ~ ' 


Y {m/2-MB-y-Sf,Sf I J,Mb-ix)Y^ 

v=±l 

1 1 


Lb — ^, 1 ,L 


ikf) 


X 


xH 


+ 


M + Ec M + E^ 


B 




M + Ec M + E_ 


B/ -i 


L,J,Lb 


[P. 


Mb j— 


V 


MB,-^l,n ^2cos^ + isin^ 


Lb — 


V2 


(25) 


where Ii,l,j,Lb. ^imLb^ PlbSIl. Plb-iIl, and involve 

radial integrals and Clebsch-Gordan coefficients. These functions are dehned 
in appendix A. The corresponding expression for the amplitude to second 
order in Ym same structure as (P5|) but the contributing terms are 

more complicated. 
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2.2 Results of the Effective Pauli Reduction 


In the previous subsection we illustrated how the effective Pauli reduction of 
the relativistic amplitude for the knock-out contribution to ( 7 ,p) reactions 
is performed to get an expansion in powers of The successive terms 

in this expansion can be reduced, in the appropriate limits, to forms that 
are equivalent to the amplitudes used in nonrelativistic calculations. This 
allows us to carry out quantitative comparisons between the relativistic and 
nonrelativistic calculations. The appropriate limits referred to above include 
i) setting the nucleon total energy equal to its rest mass ii) turning off the 
nuclear potentials in the second-order interaction Hamiltonian iii) taking 
proper account of wave function normalizations. This comparison will be 
carried out for the differential cross section as well as the photon asymmetry 
at representative energies. We shall compare the following four types of 
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calculations: 


a) Full relativistic calculations using the amplitude given by equa¬ 
tion (|lD. In the figures that follow these calculations are repre¬ 
sented by solid curves and denoted “Relativistic”. For these cal¬ 
culations the relativistic ( 7 ,p) code of Lotz has been used 


b) First order nonrelativistic calculations. These calculations are 
obtained using the amplitude of equation (pS]). Note that in 
the limit Eb and Ec —> M, the interaction Hamiltonian takes 
on a simplified form. These calculations essentially represent the 
standard nonrelativistic calculations. Comparison of these results 
with the relativistic calculations gives the essence of the differ¬ 
ence between the two approaches. These calculations are denoted 
“Pauli N.R. (First order)” in Figs. 1 and 2 and are represented 
by the dotted curves. 

c) The third type of calculation represents a nonrelativistic cal¬ 
culation carried out to second order in the inverse nucleon mass. 
This calculation then includes the interaction Hamiltonian H) , 
but with the nuclear vector and scalar potential set equal to zero 
in the interaction Hamiltonian. This interaction will then relate 
to the limit in which the photon is interacting with a free nucleon. 
Our intention here is to clarify how much improvement in the 
nonrelativistic calculations can be obtained by including second 
order effects. We shall show that the effects are not substantial. 
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These results are shown by the dot-dashed curves in Figs. 1 and 
2 and are labelled “Pauli N.R. (First -|- Second)”. We shall refer 
to these calculations in the text as ‘^medium-uncorrected' calcu¬ 
lations to signify the fact that they pertain to the limit in which 
the nuclear potentials are set to zero. 

d) The fourth type of calculation represents a nonrelativistic cal- 
culation using the full expression for Hj , i.e. with the effect 
of the nuclear medium (through the presence of the potentials) 
taken into account. These calculations are shown by the dashed 
curves in Figs. 1 and 2 and are labelled “Pauli N.R. (First -|- Full 
Second)”. They will be referred to as medium-corrected’’ non¬ 
relativistic calculations in the following text . The essence of the 
present comparison is to show the signihcance of these medium 
effects. 


Since our aim in this paper is to compare the two theoretical models, 
namely relativistic and nonrelativistic, for the knock-out contribution to 
( 7 ,p) reactions, we do not compare the resulting observables with data. We 
refer the reader to the work of Lotz and Sherif [Q for comparison of the 
results of the relativistic model with data. 

The bound state wave functions used in the calculations reported in this 
section are generated using the Dirac-Hartree potentials of Horowitz and 
Serot [^. In all the calculated cross sections given in this paper the spec¬ 
troscopic factor takes its maximum allowable value of (2 -|- 1) for both the 

relativistic and nonrelativistic calculations. 

Figure 1 shows the calculated observables for the ^®0(7,p) reaction 
with a photon energy of = 100 MeV. The residual nucleus is in a Ipi 
single hole state. The hnal state optical potential is taken from reference 
0. The cross section curves of Fig. 1(a) show large differences between 


the hrst order nonrelativistic (dotted curve) and the relativistic calculations 
(solid curve). At forward angles the nonrelativistic calculations are almost an 
order of magnitude larger than the relativistic calculations, while at backward 
angles the nonrelativistic calculations are roughly two orders of magnitude 
smaller than the relativistic calculations. Medium-uncorrected second order 
calculations improve the nonrelativistic calculations slightly only at forward 
angles (dot-dashed curve). On the other hand the medium-corrected second 
order calculations produce large changes in the nonrelativistic calculations. 
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Note the large change in the magnitnde for both small and large scattering 
angles (dashed cnrve) which brings the nonrelativistic calcnlations into close 
agreement with the resnlts of the relativistic model. 

Calcnlations of photon asymmetry shown in Fig. 1(b) also exhibit notice¬ 
able differences between the hrst order nonrelativistic and relativistic calcnla¬ 
tions. These two calculations have different shapes and magnitudes especially 
for scattering angles greater than 40°. Medium-uncorrected second order cal¬ 
culations modify the nonrelativistic calculations slightly in magnitude while 
the medium-corrected second order calculations result in noticeable changes 
in the shape and magnitude of the nonrelativistic calculations. These changes 
are such that at forward angles the nonrelativistic calculations now overlap 
the relativistic calculations and are much closer in shape and magnitude at 
backward angles. 

Figure 2 shows similar comparisons for the reaction at a higher incident 
photon energy, = 312 MeV. The hnal state global optical potentials are 
taken from Cooper et al. |^. Figure 2(a) shows that the cross section 


obtained from the hrst order nonrelativistic calculations has both different 
shape and magnitude from the results of the relativistic calculations. The 
nonrelativistic calculations lie above the relativistic calculations for angles 
smaller than 40°, whereas for other scattering angles the nonrelativistic cal¬ 
culations he below the relativistic calculations by as much as an order of 
magnitude. Medium-uncorrected second order calculations lie closer to the 
relativistic calculations only at forward angles. Medium-corrected second 
order calculations, on the other hand, are much closer to the results of the 
relativistic calculations. 

Calculations of the photon asymmetries of Fig. 2(b) show that the hrst 
order nonrelativistic results differ from the relativistic calculations in both 
shape and magnitude. Medium-uncorrected second order calculations im¬ 
prove the nonrelativistic calculations slightly at forward angles but the over¬ 
all shape stays the same as that of the hrst order. Medium-corrected second 
order calculations modify the shape and magnitude of the nonrelativistic cal¬ 
culations at all scattering angles. It might not be clear visually that these 
changes bring the nonrelativistic calculations into noticeably better agree¬ 
ment with the relativistic ones, however a chi-squared comparison does in¬ 
deed show that the medium-corrected second order calculations are closer 
to the relativistic calculations than medium-uncorrected second order calcu¬ 
lations. The important point here is that there are large diherences in the 
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calculations when medium corrections are taken into account. 

From these examples one can see that the first order nonrelativistic calcu¬ 
lations are different from relativistic calculations. Attempts to improve the 
situation through the inclusion of the second order terms in which the pres¬ 
ence of the nuclear potentials is neglected (as would normally be done in typ¬ 
ical nonrelativistic calculations) are bound to fail in bringing the results close 
to the relativistic calculations. We have seen that medium-corrected second 
order calculations (i.e. those that include the effects of the nuclear potentials 
on the interaction Hamiltonian) are much closer. We remind the reader that 
these potentials are absent in interaction terms of the usual nonrelativistic 
calculation even when higher order relativistic corrections are included [|n 


This medium modification is the important ingredient that is missing in or¬ 
dinary nonrelativistic calculations. It is this medium modification that is 
responsible for many of the differences between the two types of calculations. 
In the following section we show that the same conclusion can be reached 
through a procedure based on the Foldy Wouthuysen transformation. 


3 Foldy-Wouthuysen Transformation 

Following McVoy and Van Hove [Q many authors construct a nonrelativistic 
model of photons interacting with nuclei by performing a FW transformation 
on the relativistic Hamiltonian which involves the electromagnetic interaction 
with a free nucleon. The resulting nonrelativistic interaction Hamiltonian is 
then sandwiched between Schrodinger wave functions describing the initial 
and hnal nucleons |^. To investigate the effect of the nuclear medium on 
the FW reduction of the relativistic amplitude of the knock-out contribu¬ 
tion to (7,p) reactions, we perform a FW transformation on the relativistic 
Hamiltonian of a photon interacting with a nucleon in the presence of strong 
scalar and vector potentials. 

In the preceding discussion of the Pauli reduction scheme, our starting 
point was the distorted wave amplitude of equation (|^). A feature of this 
amplitude for practical calculation is the use of complex vector and scalar 
potentials to describe the interaction of the outgoing nucleon with the resid¬ 
ual nucleus. The unitary requirement for the FW transformation makes it 
unacceptable to work with this type of amplitude (note that in the usual 
DWBA amplitude the initial and final states of the nucleon are described by 
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different Hamiltonians). We must therefore work with a model amplitude 
in which the bound and continuum state potentials are the same and real. 
It is known that such model amplitudes are inferior in their description of 
the data in comparison to the distorted wave amplitudes. It must be noted 
however that the purpose of the present investigation is not aimed at fit¬ 
ting data; rather we are interested in features that differentiate between the 
relativistic and nonrelativistic calculations. For this purpose the restricted 
model amplitude used here is quite appropriate. 

3.1 Spin I Particle Interacting with Strong and Elec¬ 
tromagnetic Potentials 

We first summarize the FW transformation for the case in which a Dirac 
particle interacts with a general held following the procedure as given in the 
reference |^. The results will then be applied to the case in which a nucleon 
interacts with an electromagnetic held while under the inhuence of the strong 
potentials. The Dirac equation is written in the general form 

idt Ip (x) = Hip (x) , (26) 

where the relativistic Hamiltonian H can be written in terms of even and 
odd operators as 

H = pM + S + 0. (27) 

The odd operator O connects the upper component of the Dirac spinor to 
the lower component while the even operator £ can only connect either up¬ 
per or lower components. Now we perform the FW transformation on the 
relativistic Hamiltonian ( 0 ) following reference |p. After three successive 
transformations we hnd. 

(For simplicity the transformed Hamiltonian is denoted H', this is in fact the 
same as the Hamiltonian H"' in reference (p.) This result applies to any pair 
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of even and odd operators. We are interested here in a transformed Hamil¬ 
tonian containing terms which are even up to second order in the inverse of 
the nucleon mass so we will eventually drop the last term in equation (P5|). 

The relativistic Hamiltonian for a proton interacting with an external 
electromagnetic held while under the inhuence of the strong vector and scalar 
nuclear potentials can be written as 

H = a- 7r + e<h + /?[M + S{r)] + V{r) + (29) 

where tt = p — eA. The above Hamiltonian can be written in terms of the 
electric and magnetic helds E and B as 

H = a-7r-F/?[M + ^(r)]+H(r) + ^(7- E + ■ B), (30) 

where we have dehned 


S =/3cr . 


(31) 


We rewrite the Hamiltonian (^0]) in terms of even and odd operators as in 
equation (^). In this case the even and odd operators are 


£ = (3S{r) + V{r)- 


ne 

2M' 


S ■ B, 


(32) 


and 


O = a. ■ 7z -\ - 7 ■ E. (33) 

These even and odd operators are used in equation (^BD in order to get the 
FW transformed Hamiltonian for a nucleon interacting with electromagnetic 
and strong potentials. We remind the reader that the quality of the FW 
transformation depends on the assumption that the potential depths are 
small compared with the nucleon mass; V, S M. This is not the case for 
the Dirac potentials, and so we are by no means guaranteed that the FW 
transformation will yield a convergent series. We will return to this point 
below. 

The resulting Hamiltonian to second order in the inverse of the nucleon 
mass is a 4 X 4 block diagonal matrix 


H' 


\H'u 0 1 

[ 0 i/'J 


(34) 
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The upper left element of this matrix (i.e H[^) corresponds to the transformed 
Hamiltonian for positive energy solutions of the Dirac wave function. For the 
nonrelativistic limit we will use this part of the Hamiltonian which we write 
as 


H[^ = Ho + Hj, 


(36) 


where Hq involves strong potentials whereas Hj carries both strong and elec¬ 
tromagnetic interactions: we will treat the latter as a perturbation on the for¬ 
mer. These Hamiltonians are evaluated explicitly using the Coulomb gauge. 
Ho can be written to second order in 1/M as 


Ho = 


_Pl 

2M 


4M2r 

1 


+ M + S{r) + V{r) 
[S'{r)-V'{r)](7 -L 


2M2 


^(r)p2+ p[^(r)] - pT 


p^ [V{r) S{r)] 


(36) 


The interaction Hamiltonian is written in orders of 1/M as 

H, = + + 


(37) 


where for the hrst and second orders we have 


//W = 


—A ■ p 
M ^ 


2M 


(1 + k )(7 ■ ( V X A), 


= ^^(1 + 2k)|o-- V X [A]-22cr- Axpj 

+ 4^ { ■ P + 2A • p^(r) + 2S{r)c7 ■ V x [A] 

+cr. Ax V [V{r)-S{r)] 

- V{r) - S{r) ) cr • V X [A] |. (38) 


Note that in equation (^^ terms of order are dropped, and as before 
0|/]" ■ means that operator O operates only on function /. 
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3.2 Schrodinger-like Wave Functions 

The wave functions describing either the bound or continuum nucleons are 
obtained by solving the equation 


/fo^(r) = ET(r), 


(39) 


where Ho is the Hamiltonian (|3^ containing terms to second order in 1/M, 
and E is the total energy of the nucleon. Note that Hq contains a hrst 
derivative term which can be eliminated using the transformation 


T(r) = DpwW^SchW, 


(40) 


where Upw = 1 “ yv?- With this choice of Upw, the two functions ^/(r) 
and \hsch(i') have the same asymptotic form. After this transformation the 


wave equation (|39D takes the form of the Schrodinger-like equation (JT^). The 
central and spin-orbit potentials in the present case are 


Ucent (^) — 


D 


FW 


S{r) + V{r) + 


AMH 


[5'(r) + V'{r)] 


|S"(r) + \/"(r)]} 



'D' 


FWl 


f/so(r) = 


2Mr Z1pw(f) 

1 r 1 


Zlpw(r) I 4M2r 


AM 8M \D-pY>/[r) 

|S'(r) - V"(r)|' 


(41) 


The nonrelativistic amplitude for the knock-out contribution to the ( 7 ,p) 
reaction at the desired order of 1/M takes the same form as given in equa¬ 
tion (1^). Note however that the second order amplitude obtained through 
the FW scheme involves wave functions which are solutions of equation (|3^). 
As we have seen above, this equation does not have the form of the usual 
Schrodinger wave equation because it contains a hrst order derivative of the 
wave function. To be consistent with the usual nonrelativistic formalism we 
rewrite the second order amplitude in terms of the Schrodinger-like wave 
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functions \l/sch(i') introduced in equation (^0[). This requires that the inter¬ 
action Hamiltonian be modihed to 


H'j = ^ 


''FW^^^FW' 


(42) 


With this modihcation the first order terms in the interaction Hamiltonian 


will remain the same as in 
interaction Hamiltonian become 


), while the second order terms in the 


H 


( 2 )' _ 


ecu 

e 


(1 -|- 2k) {<t • V X [A ] — 2cr ■ A x V } 


4M2 I 


I 2KS{r)(T ■ V X [A] + o- ■ A X V [V{r) - ^(r)] 


- H(r) - S(r) ) (T ■ V X [A]| 


(43) 


The calculations of the amplitudes proceed in the same manner discussed in 
section 0. In particular the amplitude has the same form as that of equation 


(pTD, but now the interaction Hamiltonian is that of equation (|43|) while the 
wave functions are solutions of a Schrodinger-like equation using the central 
and spin-orbit potentials of equation (|4lD . Note that at hrst order the FW 
formalism does not produce an acceptable nucleon wave function since there 
is only a central potential and no spin-orbit potential, see equation (|4T|). The 
spin-orbit potential is borrowed from the second order and used with the hrst 
order terms to get nucleon wave functions for use in the hrst order which at 
least have a reasonable form. 


3.3 Differences Between the FW and Pauli Reduction 
Schemes 

The Pauli reduction scheme begins with the relativistic distorted wave am¬ 
plitude in which the initial and hnal nucleons are described by Dirac wave 
functions. The Dirac wave functions are solutions of the Dirac wave equation 
with diherent strong potentials for bound and continuum nucleons. 

The amplitude resulting from the Pauli reduction procedure has bound 
and continuum nucleons described by nonrelativistic wave functions. These 
wave functions are solutions of the Schrodinger-like equation (|^) with central 
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and spin-orbit potentials which in tnrn, are fnnctions of the Dirac potentials 
S{r) and V{r). 

In the FW scheme the initial bonnd and final continnnm nncleons are 
described by the same Hamiltonian. Dne to the nnrealistic potentials nsed 
to describe the continnnm nncleon, the FW calcnlations provide a toy model 
which will shed more light on the differences between the relativistic and 
nonrelativistic calcnlations. We shonld stress that nnlike the Panli redac¬ 
tion in which the wave fnnctions describing nncleons are the same in all the 
calcnlations of different orders, in the FW transformation the nncleon wave 
fnnctions are changed in each order of calcnlation dne to the contribntion of 
different terms to the wave eqnation. Differences between these two redac¬ 
tion formalisms in the free nncleon limit has been stndied by Fearing et al. 
r^ . They find that for the interaction of a real photon with a free particle 


at first order in the conpling constant the FW and Panli redactions prodnce 
the same nonrelativistic interaction Hamiltonian. 


3.3.1 Differences between the Panli and FW Hamiltonians 


In order to make instrnctive comparisons between the Panli and FW Hamil¬ 
tonians we need to simplify the Hamiltonians obtained throngh the Panli 
redaction scheme of section]^. First the potentials describing the initial and 
final nncleons are taken to be the same; in particnlar they will be the real 
Hartree potentials nsed to describe the bonnd state in the relativistic calcn¬ 
lations. In addition, the energies of the bonnd and continnnm nncleons are 
set eqnal to the nncleon mass in the Panli interaction Hamiltonians. In this 
limit the hrst order terms in the FW and Panli interaction Hamiltonians are 
exactly the same. Differences will appear in terms of order ^ and higher. 
The difference between the interaction terms of the two schemes np to second 
order, with the above modihcations of is 

(V yA)-2<T- (AxV)} 
-^{[SM + r(r)]o-. (Vx A)}. (44) 


In the special case describing the interaction of a photon with a free nncleon, 
i.e in the limit when the strong potentials S'(r) and V{r) are set eqnal to 
zero, eqnation agrees with the resnlts of Fearing et al. 


18 . Detailed 


calcnlations show that the hrst term on the right hand side of eqnation (f 
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is very small for photon energies less than a few hundred MeV, but becomes 
more important at higher energies. The second term involves the sum of 
scalar and vector potentials. At the origin the sum of these potentials is 
about -100 MeV. Note that the square of the nucleon mass (~ 1000 MeV) 
appears in the denominator of this term, so it has a very small coefficient, 
with the result that the second term on the right hand side of the above 
equation also makes only a small contribution to the transition amplitude. 
Thus the differences between the FW and Pauli interaction Hamiltonians 
appear to be small if we restrict ourselves to second order in the inverse 
nucleon mass. 

The wave functions describing nucleons in the Pauli scheme are the same 
for all orders of calculations while in the FW scheme the wave functions 
change at every order. The FW central potential at second order is similar 
in shape and energy dependence to the central potential obtained via the 
Pauli formalism. A very interesting point to note is that there is no spin- 
orbit potential at hrst order in the FW formalism: it appears at second order 
with a shape similar to that obtained in the Pauli scheme, but its magnitude 
is roughly a factor of two smaller than the spin-orbit potential arising in the 
Pauli scheme. 


3.4 Results of the FW Reduction 


In this section we compare the theoretical results of the relativistic approach 
with those of the nonrelativistic amplitudes obtained through FW reduction. 
The results are shown for three different incident photon energies. As in the 
Pauli discussion of section |2]^ we present three different types of nonrelativis¬ 
tic calculations, namely: hrst order, medium-uncorrected second order and 
medium-corrected second order. Recall that in the hrst order nonrelativistic 
calculations the wave functions are solutions of equation (|3^), ignoring all 
the second order terms in Hq (equation (|36[)) except for the spin-orbit poten¬ 
tial. In the second order calculations the wave functions are obtained from 
equation ( P5|) with all the hrst and second order termsin Hq included. The 
Hartree potentials used in all calculations (relativistic and nonrelativistic) 
are from reference |^. The graphs discussed in this section are labeled as in 
the Pauli discussion of section except for an obvious change of notation. 

Figure 3 shows the calculated observables for the {'y,p) reac¬ 
tion for a photon of energy = 100 MeV. Figure 3(a) shows the cross 
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sections. At small angles the first order nonrelativistic calculations (dotted 
curve) are about an order of magnitude lower than the relativistic calcula¬ 
tions (solid curve), while for large angles the hrst order calculations lie above 
the relativistic calculations. Medium-uncorrected second order calculations 
(dot-dashed curve) show substantial increase in the the magnitude of the 
cross sections at small scattering angles as well as some change in the shape 
of the resulting curve. Medium-corrected second order calculations (dashed 
curve) produce a noticeable change in the cross sections at backward angles. 

The photon asymmetry calculations of Fig. 3(b) also show noticeable 
differences between the first order nonrelativistic and relativistic calculations 
at backward angles. Medium-uncorrected second order calculations produce 
a change in the magnitude and the shape of the asymmetry for scattering 
angles greater than 80°. Medium-corrected second order calculations produce 
a shift towards larger angles resulting in a qualitatively similar shape to that 
of the relativistic calculations. 

Figure 4 shows the observables for the same reaction as Fig. 3 but the 
photon energy in this case is = 196 MeV. The cross section results are 
shown in Fig. 4(a), where we note that the hrst order nonrelativistic cal¬ 
culations are generally lower than the relativistic calculations by one to two 
orders of magnitude. They also fail to reproduce the dip near mid-angles. 
Second order calculations lead to a drastic change in the cross section with 
large differences due to medium corrections at both forward and backward 
angles. Medium-corrected second order calculations are in noticeably closer 
agreement with relativistic calculations compared to the medium-uncorrected 
ones. 

Similar features are observed for the photon asymmetry calculations (Fig. 
4(b)). Here again we notice large differences between the relativistic and 
hrst-order nonrelativistic calculations. Large differences also exist between 
medium-corrected and medium-uncorrected calculations. The level of agree¬ 
ment between the second order medium-corrected calculations and the rela¬ 
tivistic calculations is not the same as observed in the case of cross sections. 

The calculations shown in Fig. 5 for a photon energy of = 312 MeV 
show essentially the same qualitative features. 

One characteristic that emerges from the above discussion is that the full 
second order calculations (medium-corrected calculations) in the FW scheme 
are not as close to the relativistic results as in the Pauli reduction case at the 
same photon energy. This brings out an essential difference between the Pauli 
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and FW calculations: The wave functions in the Pauli formalism remain the 
same while different orders of the amplitude result solely from the expansion 
of the interaction Hamiltonian. The FW calculations, on the other hand, 
involve an expansion affecting both the wave functions and the interaction 
Hamiltonian simultaneously. This difference is at the root of the different 
convergence properties of the two approaches. We find that in most cases, 
by second order the Pauli expansion is quite close to the fully relativistic 
calculations, provided medium corrections are taken into account. In the 
FW scheme the level of agreement at the corresponding order is inferior, 
indicating that the convergence in this scheme is much slower than in the 
Pauli case. 


4 Conclusions 

We have described two different nonrelativistic reduction schemes of the rel¬ 
ativistic amplitude describing the knock-out contribution to ( 7 ,p) reactions. 
These reductions allow us to carry out controlled comparisons between the 
relativistic and nonrelativistic calculations of the reaction observables. In the 
Pauli formalism the relativistic S matrix is written in terms of nonrelativis¬ 
tic two-component wave functions and an effective interaction Hamiltonian. 
The effective Hamiltonian is expanded in powers of 1/ (F^ -|- M), where M is 
the nucleon mass and E is its total energy. In the limit E —> M, the hrst 
order interaction terms are exactly the same as those appearing in the usual 
nonrelativistic amplitude. The nonrelativistic wave functions in this scheme 
are solutions of the Schrodinger-like wave equation ([T^). Detailed compar¬ 
isons between the relativistic and first-order nonrelativistic predictions for 
the differential cross sections and photon asymmetries show large differences 
between the two types of calculations. The inclusion of terms to second 
order in 1/M in the interaction Hamiltonian, where medium corrections ef¬ 
fected by the nuclear potentials are left out, does not lead to any substantial 
improvement in the agreement between the relativistic and nonrelativistic 
calculations. On the other hand the expansion scheme shows explicit de¬ 
pendence in the second order terms on the nuclear potentials. When these 
medium corrections are taken into account the nonrelativistic calculations 
converge close to the relativistic results. This indicates that the essential 
difference between the relativistic and traditional nonrelativistic amplitudes. 
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is the absence in the latter of the medium modihcation of the interaction 
Hamiltonian as a consequence of the presence of the strong vector and scalar 
potentials. 

These conclusions are further supported through an analysis based on the 
Foldy-Wouthuysen transformation of the relativistic Hamiltonian describing 
a photon interacting with a nucleon embedded in the nuclear medium. The 
nonrelativistic wave functions for the bound and continuum nucleons are so¬ 
lutions of the wave equation obtained as a result of the transformation. The 
scheme leads to a nonrelativistic amplitude calculated to the desired order in 
1/M. We use these amplitudes to carry out comparisons between relativistic 
and nonrelativistic calculations in the manner described above for the Pauli 
scheme. We hud that the medium modihcations in the second order calcula¬ 
tions are important and their inclusion leads in general to better agreement 
with the relativistic calculations. However the convergence is not as efficient 
at this order as in the Pauli case. The reasons for this can be understood 
in terms of the formal differences between the structure of the nonrelativis¬ 
tic amplitude obtained using this transformation as compared to the Pauli 
reduction case. The wave functions obtained through the FW reduction are 
different at each order in 1/M, in contrast to the Pauli wave functions which 
are unchanged for all orders (recall that in the Pauli reduction only the in¬ 
teraction Hamiltonian is expanded). 

The basic result of the present work is that standard nonrelativistic cal¬ 
culations of the knock-out amplitude do not properly take into account the 
strong medium modihcations of the interaction Hamiltonian. We have clari- 
hed this point through a comparison based on nonrelativistic reduction of the 
relativistic amplitude using both the Pauli and Foldy-Wouthuysen reduction 
schemes. 
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Appendix A 


In this appendix we give the explicit form of the six functions involving 
Clebsch-Gordan coefficients and radial integrals which are introduced in the 


amplitude of equation (|25|). 

= (Ls,/;0,0 I L,0) [ r'^drfB{r)ji{k^r)fLj{r), 


H,L,J,Le 


^Lb+ 1,1,L 


V 






^L,J,Lb 


<5,Lb = {L,l/2-MB + fi,fi\J,MB + 2fi) 

X {Lb, 1/2; Mb + p, —p | Jb, Mb), 


777 -~r(-h_B + 1 ; 1 ; Mb fi — i', i' I Lb, Mb + p-) 

2Lb + i 

{Lb 1,1 Mb + p — z/, 0 | L, Mb + p — i^){L, /; 0, 0 

{Lb — 1,1; Mb + p — z/, z/ | Lb, Mb + p) 


Y 2Lb + 1 

X {Lb — 1, / : Mb + p — z/, 0 | L, Mb + P ~ z/) 
x(Z/,/;0,0 I Lb — IjO) 

X J fLj{r)ji{k^r) r^dr, 

(L, 1/2; Ms - p - p p I J^Mb-v) 

X {Lb, 1/2; Mb — p, p | Js, Ms). 


^i,l,j,Lb ~ {Lb, 1 / 2 ; Mb + fi,-fi \ Jb, Mb) 

X {L, 1/2; Mb “1“ P; P | '.f; Mb ~\~ 2p) 
X (L, /; Mb + p, 0 | Ls, Mb + p), 


+ 1 , 0 ) 
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Figure Captions 


Fig. 1: Differential cross section (a), and photon asymmetry (b) for the re¬ 
action ^®0(7, pY^N at = 100 MeV. 

Solid curve - full relativistic calculations. 

Dotted curve - nonrelativistic calculations using the hrst order Hamiltonian 
hY^ of equation (plQl). 

Dot-dashed curve - second order nonrelativistic calculations (neglecting the 
nuclear potentials in from equation (|I^)). These are referred to as 
medium-uncorrected second order calculations in the text. 

Dashed curve - second order nonrelativistic calculations using the full Hj . 
These are referred to as medium-corrected second order calculations. 


Fig. 2: Differential cross section (a), and photon asymmetry (b) for the 
reaction of Fig. 1 but with = 312 MeV. Curves labelled as in Fig. 1. 


Fig. 3: Differential cross section (a), and photon asymmetry (b) for the 
reaction ^®0(7, pY^N at = 100 MeV. 

Solid curve - full relativistic calculations. 

Dotted curve - nonrelativistic calculations using the hrst order FW Hamilto¬ 
nian hY^ of equation (|5BD and nucleon wave functions obtained from equation 
(|39D using the hrst order and spin-orbit terms of equation (|36D . 

Dot-dashed curve - second order nonrelativistic FW calculations (neglecting 
the nuclear potentials in iF) from equation (|i3|)), and full second order 
potentials to generate the wave functions. These are referred to as medium- 
uncorrected second order calculations in the text. 

Dashed curve - second order nonrelativistic FW calculations using the full 
Hj . These are referred to as medium-corrected second order calculations. 


Fig. 4: Diherential cross section (a), and photon asymmetry (b) for the 
reaction of Fig 3 but with E^ = 196 MeV. Curves labelled as in Fig. 3. 

Fig. 5: Diherential cross section (a), and photon asymmetry (b) for the 
reaction of Fig 3 but with E^ = 312 MeV. Curves labelled as in Fig. 3. 
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